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<^ , Abstract. I present a review by epoch of baryons in the intergalac- 

■ tic medium (IGM), from the first star until today. Recent observations 

indicate a protracted period of reionization, suggesting multiple popula- 
tions of reionizers; detection of these z ^ 6 sources is a key goal that is 



(N 
> 
00 



CO 



now coming within reach. The optical Lyman alpha forest (2 <J z <^ 4) 
■^j" \ is well-described by the Fluctuating Gunn-Peterson Approximation, but 

puzzling observations of associated metal lines and nearby galaxies may 
provide insights into galactic feedback processes. IGM studies at z ^ 1.5 
are progressing rapidly thanks to ultraviolet absorption line studies of 
both Lya absorbers and the warm-hot intergalactic medium. Cosmo- 
lO ■ logical hydrodynamic simulations have played an integral part in these 

advances, and have helped to shape our understanding of the IGM at 
each epoch. Working jointly, observations and theory continue to expand 
our knowledge of the IGM as the earliest stages of galaxy formation and 
the dominant reservoir of baryons at all redshifts. 
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• i— i . 1. Introduction 

X 

H . ... 

All baryons in the Universe were originally in the intergalactic medium (IGM). 

Today, the majority of baryons remain there. Understanding the evolution of 
the IGM is therefore a critical component for a complete theory of how baryons 
evolve in the Universe, from the smooth Dark Ages to present-day galaxies and 
large-scale structure. In this review I will highlight some recent work towards 
understanding the evolution of the IGM's various components, their physical 
state and connection to galaxies, with particular emphasis on the contributions 
and insights provided by numerical simulations. 

This review is divided into four sections. §2 presents a simulation-based 
overview of the various baryonic components in the Universe. §3 discusses some 
of the latest theoretical and observational results probing the epoch of reioniza- 
tion. §4 presents a short history of optical Lya forest studies, and expounds on 
two of its more puzzling aspects: The presence of metals in the diffuse IGM, 
and the effects of galaxies on their surrounding gas. §5 highlights recent de- 
velopments in probing the z 1.5 IGM, including prospects for detecting the 
so-called missing baryons. §6 provides a short summary. 
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Figure 1. Left panel: The cosmic phase diagram at z = from a 
ACDM simulation by Springel & Hernquist (2003). Phases described 
in text are demarcated. Right panel: The evolution of the mass fraction 
of baryons in each phase from z = 7 — ► 0. 



2. The History of Baryons in the Universe 



Baryons in the Universe can be divided into four phases, according to their 
density and temperature (e.g. Dave et al. 1999): (1) "Condensed": Stars 
and cold gas in galaxies. This is the most readily detectable baryonic phase. 
(2) "Hot": Gas bound in clusters. These baryons are seen via their X-ray 
emission, which generally requires that they have T £ 10 7 K. (3) "Diffuse": 
Low-density gas whose dynamics are governed by the underlying dark matter 
distribution. This phase gives rise to the majority of Lya absorbers seen in 
distant quasar spectra. (4) "Warm-Hot": Mildly shock-heated, uncollapsed gas 
with 10 5 £ T £ 10 7 K. This gas is difficult to detect in absorption due to its 
high ionization state, and difficult to detect in emission due to its low density. 
Hence these are the "missing baryons" that have drawn considerable attention, 
as we will discuss in §5. 

Figure 1 (left panel) shows the z = cosmological phase diagram, i.e. 
a density-temperature plot of gas in a representative volume of the Universe, 
drawn from a TreeSPH simulation by Springel & Hernquist (2003, their 100 Mpc/h 
G4 run). They use kinetic feedback and a multi-phase ISM model to obtain a 
converged star formation history that reproduces the observed stellar density at 
z = 0. Divisions between the four phases above are indicated. At the lowest 
densities, the diffuse phase shows a tight density-temperature relation (discussed 
in §4.2). Moving to mildly overdense regions allows gas to compress and shock 
heat, producing a scatter in temperature that grows with density through the 
warm-hot phase. As gas falls into the largest bound systems it virializes at the 
halo temperature and enters the hot phase. At larger densities, the cooling time 
becomes short and gas is able to cool into the condensed phase and form stars. 
The upturn at p/p > 10 6 arises from the details of feedback implementation in 
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this simulation, and the floor at ~ 10 4 K in the condensed region arises because 
no molecular hydrogen or metal line cooling is included in this run. 

The right panel of Figure 1 shows the mass fraction in each phase from 
z ~ 7 — > 0. At z ^ 3, virtually all baryons are in the diffuse phase. As 
structure forms, intergalactic gas in the filaments surrounding galaxies is able 
to shock-heat, and the warm-hot phase begins to grow. The condensed phase 
increases rapidly up to z ~ 2, then the growth slows considerably, following 
the observed cosmic star formation history. Cluster formation governs the hot 
phase; this becomes significant only at late times when large potential wells are 
able to form. By the present day, ~ 45% of the baryons reside in the warm-hot 
phase. This conclusion is broadly consistent with that from a number of other 
cosmological hydrodynamic simulations, using a variety of numerical techniques, 
resolutions, and input microphysics (Dave et al. 2001). Thus these predictions 
for the evolution of baryons are generic to CDM models, and are driven primarily 
by the gravitational growth of structure (Valageas h Silk 2001). 

3. The Reionization Epoch 
3.1. The First Stars 

The IGM came into being, by definition, when the first galaxies formed. At 
this time, the Universe was essentially metal- and dust-free, with (presumably) 
negligible magnetic fields. Thus the first star formed in a relatively simple envi- 
ronment, from the gravitational collapse of a primordial gas cloud (Tegmark et 
al. 1997). Ultra- high resolution simulations by Abel, Bryan & Norman (2002) 
and Bromm, Coppi & Larsen (2002) show that a star produced in this environ- 
ment is likely to be a VMS (very massive star, M ^ 100M Q ), since the lack of 
H2 and metal cooling agents results in a large Jeans mass. 

The calculation of the second star is significantly more complex. Due to the 
exceptionally biased nature of the first collapsed object (an ~ 8a perturbation), 
the second star would be expected to form close to the first star. However, the 
first star emits a copious amount of ultraviolet radiation during its <J 3 million 
year lifetime, which can sufficiently heat surrounding gas so as to inhibit subse- 
quent collapse. If the collapse is inhibited beyond the lifetime of the first star, 
then the resulting supernova may have a dramatic impact on the second star's 
gas cloud. VMS's were originally believed to collapse directly to black holes 
with no ejectae; however, calculations by Heger & Woosley (2002) show that 
in a specific stellar mass range (« 140 — 260M Q ), the resulting pair instability 
supernova can expel a large fraction of its material into the surrounding IGM. 
Not only would the mechanical wind energy evacuate surrounding gas, but the 
accompanying metals would pollute the IGM. If the pollution reaches a modest 
value of ~ 10~ 4 Z Q , then subsequent star formation from that gas will not be 
VMS's but will instead follow a Population II IMF (Schneider et al. 2003). 

From these considerations, it might appear that the first few massive stars 
would serve to quickly pollute the IGM and transition to a lower-mass IMF 
incapable of reionizing the Universe until late times. This cannot be the whole 
story, though, since WMAP observations of the temperature-polarization cross 
correlation indicate a Thomson optical depth of r e = 0.17±0.04, implying a high 
reionization redshift of « 11 — 30 (95% confidence; Kogut et al. 2003). While 
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the error bars are still large (and will improve with accumulating WMAP data) , 
if taken at face value they suggest a sizeable population of efficiently-emitting 
VMS's in place at early times (Cen 2003b). 

3.2. The Minihalo Epoch 

The epoch of minihalos, i.e. bound systems with virial temperatures below 
10 4 K that can only cool via molecular hydrogen line emission in a primordial 
gas, is a crucial one for the reionization history of the Universe. Minihalos 
are so ubiquitous at z ~ 10 — 15 that they are expected to have a covering 
fraction of unity around all halos having T v i r i a \ > 10 4 K that are likely to be 
ionizing sources, thereby providing a large sink for reionizing photons (Shapiro, 
Iliev & Raga 2003). However, if minihalos are able to form stars themselves, 
they may provide a large source of reionizing photons instead. H2 production in 
primordial gas is insufficient to rapidly cool minihalos, but the presence of metals 
or X-rays that can stimulate H2 production could alter this scenario (Haiman, 
Abel & Rees 2000; though see Venkatesan, Giroux & Shull 2001). An early 
population of black holes (Madau et al. 2003) may be a candidate for producing 
such hard X-rays. However, Oh & Haiman (2003) argue that even if a minihalo 
forms an initial burst of stars, the resulting increase in entropy in the relic 
H II region would prevent further gas collapse; thus they argue that minihalos 
cannot contribute significantly to r e measured by WMAP. Given the conflicting 
theoretical predictions, direct observations of minihalos are highly desireable. 
Prospects include seeing them in 21cm emission with LoFAR or ALMA (Iliev 
et al. 2002, 2003), or in 21cm absorption to background sources (Furlanetto & 
Loeb 2003) if such sources exist. 

Since the Universe at z ^ 10 is sufficiently dense to recombine after being 
initially ionized, and feedback from the first generation of stars may suppress 
subsequent star formation, it is possible that the Universe underwent two epochs 
of reionization (Cen 2003; Wyithe & Loeb 2003a). The first would occur at 
z £ 15 from VMS's (Oh et al. 2001; Wyithe & Loeb 2003b), and the second 
around z ~ 6 from Population II stars in protogalaxies. This would reconcile 
the WMAP results with observations of a sharp rise in the IGM neutral fraction 
seen in spectra of z ;> 5.7 Sloan Digital Sky Survey quasars (White et al. 2003). 
Sloan data also shows that the numbers of z ~ 6 quasars are too small to 
provide the necessary photons to perform this final stage of reionization (Fan et 
al. 2003), thus star forming regions in protogalaxies are the likely culprit. An 
observational signature of multiple reionizations may lie in 21cm emission from 
the IGM (Furlanetto, Sokasian & Hernquist 2003). 

3.3. Detecting the Reionizers at z ^ 6 

Observing the sources responsible for the final reionization of the Universe is 
one science driver for the James Webb Space Telescope (JWST), scheduled to 
lauch in 2011. Broad-band color selection can be used to identify candidates, 
but detailed spectroscopic follow-up will be challenging from space, and may 
require novel approaches in ground-based infrared spectroscopy. 

Perhaps surprisingly, it may be more efficient to search for z ^ 6 star forma- 
tion using ground-based narrow band searches for Lya emission. Cowie (priv. 
comm.) has confirmed several z 6.6 sources this way, so Lya emission is in fact 
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Figure 2. Left panel: 32-hour "observation" with 8m telescope of 
Lya emission from z = 8.2 galaxies in a ACDM simulation. Square 
inset shows a close-up of dashed box, the most heavily populated re- 
gion. Right panel: Lya luminosity function for z = 8.2 galaxies in two 
IMF// esc scenarios , with detection limits for 8m and 30m telescopes 
with two narrow-band filters. Figures from Barton et al. (2003). 



detectable prior to the end of reionization. Barton et al. (2003) used a hydro sim- 
ulation to show that with optimistic (but not outrageous) assumptions regarding 
the IMF and escape fractions, it may be possible to detect pre-reionization ob- 
jects with existing 8-10m telescopes, using an R ;> 100 narrow-band filter in 
a night sky window corresponding to Lya emission at z = 8.2. Figure 2 (left 
panel) shows a simulated 32-hour observation with a VMS-based IMF and fa- 
vorable escape parameters from Santos (2003). These include having sources 
blow a local ionization bubble (Haiman 2003), having the emission be centered 
redwards of the protogalaxy's redshift, and no dust attenuation. Under these 
(admittedly fortuitous) conditions, a fair number of sources would be detectable. 

While detection with an 8m telescope will require some cooperation from 
nature, detection with a 20m+ class telescope would be straightforward even for 
fairly conventional IMFs and escape parameters. The right panel of Figure 2 
shows the Lya luminosity function predicted from the Q5 simulation of Springel 
& Hernquist (2003), having 2 x 324 3 particles and a 10 Mpc/h (comoving) box 
size. For a mildly top-heavy IMF as might be expected for Population II stars, 
the Lya luminosities would be ~ 6x lower than the "plausible" model. This 
would preclude detection with an 8m class telescope, at least without a large- 
area survey, but could be seen with a 20-30m class telescope. Such a telescope 
would also be required for spectroscopic follow-up of JWST broad-band selected 
high redshift galaxies, which is crucial for constraining their physical properties. 
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4. The Optical Lya Forest 

4.1. History of Cloud Models 

Redshifted intervening weak neutral hydrogen absorption, i.e. the "Lya forest" , 
was first detected in quasar spectra in the 1960's (e.g. Schmidt 1965; Lynds & 
Stockton 1966). Their cosmological origin was clearly demonstrated by Sargent 
et al. (1980), who also found a substantial number of associated C IV absorbers. 
Models for absorption systems trace back to Gunn & Peterson (1965) and Bah- 
call & Salpeter (1965), who forwarded the idea of a smooth, highly ionized IGM 
causing a uniform "Gunn- Peterson trough" (as it was later called), punctuated 
by discrete clouds of neutral hydrogen. This model persisted, in various forms, 
into the 1990's, with individual cloud confinement provided by ambient pressure 
from a hot IGM (Sargent et al. 1980; Ikeuchi & Ostriker 1986), gravity (Melott 
1980; Black 1981), CDM minihalos (Rees 1986; Ikeuchi 1986), and combinations 
thereof (Petitjean et al. 1993). With ever-improving observations, it became 
clear that no model could fit Lya data without a great deal of uncomfortable 
tuning. The advent of Keek's HIRES spectrograph, together with new insights 
from hydro simulations, paved the way for the final demise of "cloud" models. 

4.2. The Fluctuating Gunn-Peterson Model 

In the early 1990's, an alternative origin for Lya absorbers was being formulated 
from rapidly advancing cosmological hydrodynamic simulations. It was noted 
that the Cosmic Web must necessarily contain gas ionized by a metagalactic 
flux from quasars. The resulting H I absorption along any line of sight had 
statistical properties remarkably similar to that of observed Lya systems (e.g. 
Cen et al. 1994; Zhang, Anninos & Norman 1995). In particular, the column 
density distribution (i.e. the number of systems per unit column density per 
unit redshift) from such models naturally reproduced the observed power law 
over ~ 10 orders of magnitude in Nm (Hernquist et al. 1996; Dave et al. 1997). 

From these simulations a new paradigm arose for the origin of the Lya for- 
est. Instead of arising in individual clouds, Lya absorption is analogous to a 
Gunn-Peterson trough that fluctuates as it passes through over- and under-dense 
regions of the Universe. The absorbing gas is expected to have a tight density- 
temperature relation set by the balance between adiabatic cooling due to Hubble 
expansion and photoionization heating by the metagalactic UV background (Hui 
Sz Gnedin 1997). This "equation of state" of the IGM can be constrained obser- 
vationally (Schaye et al. 2000; Choudhury, Padmanabhan, k, Srianand 2001). 
Thus a given density corresponds to a given temperature, which corresponds 
to a given ionization fraction, which in turn corresponds to a given H I optical 
depth. The simulation paradigm can thus be encapsulated in the "Fluctuating 
Gunn-Peterson Approximation" (FGPA; Croft et al. 1998), in which thi oc p 1 ' 6 . 
The density p is that of either the dark matter or the gas, since they trace each 
other quite well in these low-density, gravitationally-dominated systems. The 
constant of proportionality is inversely related to the H I photoionization rate, 
and also depends on the IGM temperature and cosmological parameters (see 
Croft et al. for precise form). 

The FGPA has proven to be a valuable tool for interpreting Lya forest data. 
For instance, according to the FGPA, the mean baryon density il^ can be simply 
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related to the mean optical depth of H I absorption. Rauch et al. (1997) mea- 
sured thi, and together with a minimal estimate of the H I photoionization rate 
from quasars, determined Qj, ^ 0.017/i -2 . At the time, this determination was 
considered quite high, though it agreed with initial results from [D/H] studies 
(Tytler, Fan & Buries 1996). WMAP results (Spergel et al. 2003) have now 
confirmed the Rauch et al. value. Another application of the FGPA is to inter- 
pret metal lines observations (discussed in §4.3), using the conversion between 
thi and (p, T) to obtain accurate, pixel-by-pixel ionization corrections. 

A more ambitious application of the FGPA, and the one for which it was 
originally proposed, is the reconstruction of the matter power spectrum (Croft 
et al. 1998). Since the optical depth is a one-dimensional, nonlinear map of 
the underlying dark matter distribution, it can be inverted to yield the matter 
fluctuation spectrum. The cosmological parameters derived from the latest such 
analysis at z = 2.7 (Croft et al. 2002) are consistent with the WMAP results at 
z = 1089 (Spergel et al. 2003) and large-scale structure studies at z ^ 0.5, but 
uniquely probes an intermediate redshift regime. WMAP and the Croft et al. 
data together show a mild preference for a running tilt of the primordial power 
spectrum, but this is controversial (Seljak, McDonald & Makarov 2003). 

4.3. Metals in the Diffuse IGM 

High resolution spectroscopy on 8- 10m class telescopes has opened up a new win- 
dow for studying metal absorption in low-density regions of the Universe. Using 
Keck/HIRES, Songaila h Cowie (1996) found that most Lya systems at z ~ 3 
with 7V"hi > 10 15 cm -2 , and half the systems with Nm > 10 145 cm _2 , showed 
associated C IV absorption. This is fairly surprising, since the latter column 
density limit corresponds to matter overdensities of a few at those redshifts, 
which presumably lie well outside of galaxies. Thus some mechanism must be 
invoked to transport metals across cosmological distances. Constraints on such 
mechanisms can be provided by examining the density dependence of the metal- 
licity, Z{p): A strong gradient would imply late, local enrichment, while a more 
uniform distribution would suggest an early epoch of enrichment. 

There are two approaches to studying Z(p): One is to trace a single ion 
over a range of IGM densities, and the other is to examine different ions that 
trace different densities. The idea behind the latter is shown in the left panel of 
Figure 3: The curves show the Line Observability Index (LOX; Hellsten et al. 
1998) for C IV and O VI, as a function of N m (which, via the FGPA, yields 
the ionization conditions). The curves shown use a metallicity [C/H]= —2.5 
and a halo star abundance pattern with [0/C]= 0.5. The LOX translates to an 
equivalent width, in mA, of a "typical" associated metal absorber. The solid 
lines show the LOX for a Haardt & Madau (1996) background, while the dashed 
lines show the effects of attenuating this background by xlO above 4 Ryd in 
order to mimic the ionization conditions prior to He II reionization (e.g. Heap 
et al. 2000). Notice that a comparison of these two ions can provide not only an 
estimate of Z(p), but also constrain the shape of the ionizing flux. The horizontal 
lines denote approximate detection limits for Keck/HIRES, for within the Lya 
forest (dashed; appropriate for O VI) and redwards of it (solid; for C IV). From 
this, it is predicted that C IV would routinely be seen for Nm > 10 14 ' 5 cm~ 2 , 
as observed, but at lower densities it is increasingly ionized to C V. Meanwhile, 
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Figure 3. Left panel: Line Observability Index (LOX) for C IV and 
O VI at z = 3, for a quasar background (solid) and a pre-He II reioniza- 
tion ("cut") background (dashed). Right panel: Apparent optical depth 
ratios for pixels in Q1422 spectrum with significant C IV absorption. 
Comparison with simulated spectra clearly indicates the presence of 
associated O VI absorption. Figures from Dave et al. (1998). 

the lowest density gas is probed nicely with O VI absorption. Systems with 
Nhi ps 10 14 cm~ 2 correspond to gas at the mean density of the Universe at 
z = 3, so in principle O VI absorption can even probe into voids. Note that 
O VI here is photoionized, not collisionally ionized, consistent with observations 
of O VI absorbers at z £ 2 (Carswell et al. 2002). 

The practical difficulty with using O VI is that its doublet lies at rest 
wavelengths 1032&1037A, firmly within the dense Lya forest, so it is quite dif- 
ficult to identify and confirm the doublet lines. Instead, Dave et al. (1998) 
followed Songaila (1998) in taking a pixel-based approach to identify an excess 
of absorption over the expected contribution from H I at the wavelength of as- 
sociated O VI. The so-called optical depth ratio technique yields the plot shown 
in the right panel of Figure 3: The histogram of apparent optical depth ratios 
Novl/Nciv from a Keck/HIRES spectrum of Q1422+2309 is compared to pre- 
dictions from a hydro simulation with and without oxygen. This shows that O VI 
is clearly present in C IV systems, at a level consistent with a quasar-dominated, 
post-He II recombination ionizing background and an alpha-enhanced abun- 
dance pattern. Dave et al. also found no O VI in weaker H I systems that 
showed no C IV. From these results, Dave et al. reached two main conclusions: 

1. A comparison of O VI in systems with and without C IV indicates a 
metallicity gradient of dZ/dlogp ^ 0.5. 

2. The presence of O VI means that He II was reionized over at least half the 
redshift path length covered by Q1422's Lya forest from z = 3.5 — > 3.1. 

Though much larger data sets now exist, this work has not been followed up due 
to the lack of availability of Keck/HIRES quasar spectra. 
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Meanwhile, VLT's UVES spectrograph has now obtained large samples of 
high-redshift quasar spectra. Aguirre et al. (2002) developed a somewhat differ- 
ent pixel-based analysis technique called the Pixel Optical Depth (POD) method 
to optimally extract metallicity information in low-density regions. Using 9 
UVES spectra and hydro simulations for calibration, Schaye et al. (2003) quan- 
tified the variation in IGM metallicity with both density and redshift, employing 
the simpler, more robust idea of measuring C IV absorption over a sizeable range 
of H I optical depths. They obtained [C/H]= -3.4718$ + 0.08ljj;?jj(z - 3) + 
0.65to:i°(log(& - 0.5), over -0.5 £ log 5 £ 1 and 2 £ z £ 4. Thus they find a 
modest gradient with density (consistent with Dave et al. 1998), but little evo- 
lution of the IGM metallicity. A similar lack of evolution in C IV was also found 
by Songaila (2001) for z = 5.5 — > 2. Aguirre et al. (2003) extended this method 
to Si IV, which generally traces higher densities than C IV, and confirmed the 
lack of evolution in a new ion, in addition to finding intriguing evidence for large 
alpha enhancements ([Si/C]« 0.77) in the IGM. These results generally point 
towards an early epoch of enrichment at z >> 3. 

Conversely, recent observations of superwinds from Lyman break [z ~ 3) 
galaxies (Pettini et al. 2002) argues in favor of a local enrichment mechanism. 
Adelberger et al. (2003) find a curious lack of H I absorption in the immediate 
vicinity of galaxies at impact parameters <J 0.5 comoving Mpc/h (see the solid 
points in the left panel of Figure 4), which they interpret as further evidence for 
strong winds. Adelberger (these proceedings) also found a strong enrichment 
of C IV along lines of sight passing near galaxies, providing direct evidence for 
in situ enrichment. Thus the exact manner and epoch at which the IGM was 
enriched remains uncertain. 

4.4. Absorption Near Galaxies at z ~ 3 

The galaxy-IGM connection explored by Adelberger et al. (2003) also presents 
a mystery. The dramatic drop in absorption near galaxies is opposite to what is 
expected from hierarchical simulations, since matter clustering tends to increase 
absorption close to galaxies (Kollmeier et al. 2003). The left panel of Figure 4 
shows the Adelberger et al. data (labelled "ASSP") compared to a hydro sim- 
ulation having no winds/photoionization (triangles), and various ad hoc model 
with complete evacuation of H I out to R w = 0.75 (squares), 1 (pentagons), and 
1.5 (hexagons) comoving Mpc/h, and another model where the evacuation vol- 
ume is proportional to the galaxies' baryonic mass (stars). Only in the extreme 
case of complete evacuation out to 1.5 Mpc/h does the model come within the 
error bars of the data. The radius would have to be even larger if the evacuation 
was (more realistically) non-spherical or incomplete. Large evacuation radii are 
required because peculiar velocities can shift absorption from large distances to 
the galaxies' redshifts (Desjacques et al. 2003). Models based on local photoion- 
ization by star formation fail by a considerable margin to produce the required 
evacuations, so winds are the logical explanation, but the required energies are 
quite large. Similar conclusions were found by Croft et al. (2002). 

The fact that only three quasar-galaxy pairs contribute to Adelberger et al.'s 
innermost point raises the possibility that the lack of absorption is a statistical 
fluke. To examine this, Kollmeier et al. sampled trios of galaxy-absorber pairs 
in various models to determine the frequency of obtaining a flux decrement as 
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Figure 4. Left "panel: Lya flux decrement as a function of radius from 
galaxies at z ~ 3 for various evacuation models described in text, com- 
pared with Adelberger et al. data. Right panel: Statistical probability 
of reproducing the innermost Adelberger et al. point, in various model 
scenarios. Figures from Kollmeier et al. (2003). 



low as observed. Figure 4, right panel, shows that it is exceptionally rare for 
our fiducial simulation, even when a random redshift error of a v = 300 km/s 
is included. Even for our most extreme evacuation model there is barely a 5% 
chance of reproducing the data. 

It remains unclear what physics is required to reconcile the models with 
this data. Kollmeier et al. considered photoionization from AGNs, "filling in " 
of Lya absorption from cooling radiation (Fardal et al. 2001), non-equilibrium 
effects, and other ideas, but none appeared promising. Obviously a larger sample 
of these systems would be helpful, and has already been obtained at z ~ 2 
(Steidel, priv. comm.). If these observational results are borne out, it may 
require a radical change in our understanding of feedback processes. 



5. The Low- Redshift Intergalactic Medium 
5.1. The Lya Forest at 1.5 

The launch of HST with its ultraviolet capabilities heralded a new age in IGM 
studies, as the Lya absorption line (and various UV metal features) became 
routinely observable for systems with z <^ 1.5. The Quasar Absorption Line 
Key Project compiled over 80 spectra (Jannuzi et al. 1998), and quantified the 
statistics of low-z absorbers with great precision (Weymann et al. 1998), even 
though the Faint Object Spectrograph (FOS) data could not resolve individual 
Lya absorbers. A major surprise was that any absorbers were found at all, 
since an extrapolation from the rapidly-declining population at z ~ 4 — > 2 
would predict virtually none (Morris et al. 1991). 

Dave et al. (1999) extended hydrodynamic simulations to z = using the 
first distributed-memory parallel TreeSPH code (Dave, Dubinksi, & Hernquist 
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Figure 5. Left panel: Column density distribution of Lya absorbers 
at z = 0.17 from PG0953+415 and H1821+643 (solid), vs. ACDM 
simulation (dashed). Turnover is due to noise; higher-S/N simulated 
spectra (dotted) show no break in power law. Right panel: Linewidth 
(^-parameter) distributions. Thermal contribution to linewidths (dot- 
ted) is significant at low-z. Figures from Dave & Tripp (2001). 

1997), and found that the statistical properties of simulated absorbers were in 
good agreement with Key Project results. Furthermore, they argued against 
a second population of slowly-evolving absorbers dominating at low-z (Bahcall 
et al. 1996), and instead claimed that the change in evolution of the absorber 
population was caused by a drop in the amplitude of the quasar-dominated 
metagalactic flux, which increased the neutral fraction and balanced the drop 
in gas column density due to cosmic expansion. In other words, low-z Lya 
absorbers, at least the weaker ones, are physically analogous to their high-z 
counterparts, albeit shifted to somewhat lower column densities. 

The Space Telescope Imaging Spectrograph (STIS), installed in 1997, pro- 
vided the first fully resolved glimpse at the local Lya forest. Careful comparisons 
versus hydro simulations continued to show excellent agreement with data (Dave 
&; Tripp 2001; see Figure 5). This suggests that with a large sample of weak 
absorbers, as should be provided by the Cosmic Origins Spectrograph beginning 
in 2005, many of the analyses done at high-z could be performed at z <; 1.5. 
Hence COS offers great promise for constraining the metagalactic ionizing flux 
(Dave & Tripp 2001), the IGM metallicity, and possibly even the matter power 
spectrum at z ^ 1 on mildly nonlinear scales. 

5.2. The Warm-Hot Intergalactic Medium 

Despite the physical similarities in Lya absorbers, the low-z IGM is considerably 
more complex than the high-z one. As shown in Figure 1, at z > 3 the diffuse 
IGM contains most of the baryons in the Unvierse, but this is not true at low 
redshifts. In particular, the elusive Warm-Hot Intergalactic Medium (WHIM) 
is predicted to contain a third to half of all baryons today (Dave et al. 2001). 
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Figure 6. Left panel: The Warm-Hot Intergalactic Medium in a simu- 
lation by Cen & Ostriker. Green through red shows T gas = 10 5 — ► 10 7 K. 
Figure from Dave et al. (2001). Right panel: Line density of O VI ab- 
sorbers seen in a cosmological hydro simulation (lines) with various 
assumption regarding how metals trace density. Data points show ob- 
servations of Tripp & collaborators. Figure from Chen et al. (2003). 



The WHIM has drawn widespread interest recently because it may be the 
repository of the so-called "missing baryons". Fukugita, Hogan & Peebles (1998) 
performed an inventory of baryons observed today and found roughly a half to 
be missing, compared with high-z Lya forest measurements or the now well- 
determined Oft. Cen & Ostriker (1999) suggested that they reside in a warm-hot 
state, with 10 5 < T < 10 7 K, making them difficult to detect. Dave et al. (2001) 
examined the physical properties of WHIM gas in more detail, finding that it 
generally resides in filamentary large-scale structures, well outside of virialized 
objects (Figure 6, left panel). The existence of such gas in large quantities 
was at first surprising, since if it is all placed in virialized halos then the soft 
X-ray background would be overproduced by orders of magnitude (Pen 1999; 
Wu, Fabian Sz Nulsen 1999). Instead, the low densities predicted by simulations 
result in a much weaker background, broadly consistent with observational limits 
(Croft et al. 2001; Phillips, Ostriker & Cen 2001). 

Detecting the WHIM in absorption is challenging, since H I is a poor tracer 
of gas at these temperatures. However, O VI is promising because it has a 
collisional ionization maximum at ps 3 x 10 5 K. Tripp, Savage & Jenkins (2000) 
detected four intervening O VI systems at z <^ 0.3, and obtained a lower limit 
on the baryon density contained in such systems to be ;> O.lfJf,. Because of 
the uncertainties in ionization corrections, metallicities, and gas densities, it is 
perhaps more straightforward to make predictions from simulations regarding 
the expected frequency of O VI absorbers. The right panel of Figure 6 shows 
such a comparison from Chen et al. (2003). Observations (solid points) are 
reproduced nicely with a metallicity of roughly one-tenth solar for WHIM gas, 
a metallicity that is somewhat higher than Schaye et al. (see §4.3) at WHIM 
densities (5 ~ 30), extrapolated to z ~ 0. Cen et al. (2001) and Fang & 
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Bryan (2002) obtain similar results. O VI absorption has also been detected in 
association with the Local Group (Nicastro et al. 2003). 

It is also possible to detect WHIM absorption in the X-ray band, par- 
ticularly using the O VII and OVIII Lines (Chen et al. 2003; Fang, Bryan 
& Canizares 2003). Mathur, Weinberg & Chen (2003) claim to have detected 
WHIM O VII absorption in a Chandra spectrum of H1821, though XMM follow- 
up did not unambiguously confirm it. 

Detecting WHIM emission, while even more difficult, may ultimately prove 
the best approach for obtaining a complete map of missing baryons in the local 
Universe. Zappacosta et al. (2002) found a diffuse soft X-ray structure associ- 
ated with a galaxy overdensity at z = 0.45, but a more systematic search will 
have to await new satellites. The proposed Missing Baryon Explorer (MBE) 
would detect a significant number of WHIM sources in six X-ray emission lines 
(Fang et al. 2003). The Japanese have a similar proposal called Diffuse Inter- 
galactic Oxygen Surveyor (DIOS). Emission in lower ionic species such as O VI 
and C V would detect the bulk of WHIM gas at lower temperatures, and the pro- 
posed Spectroscopy and Photometry of the IGM's Diffuse Radiation (SPIDR) 
satellite was aiming to do so until it was cancelled in the design phase. 

6. Summary 

Observations of the IGM have progressed rapidly in recent years, which in turn 
has spurred dramatic advances in our theoretical understanding of the various 
baryonic components in the Universe. Still, many questions crucial to our un- 
derstanding of galaxy and structure formation remain unanswered, including: 

• What are the sources that reionize the IGM, and how can we detect them? 

• How and when did metals get into the IGM? 

• What effects do galactic feedback have on surrounding IGM gas, and how 
does it change with redshift? 

• How do we obtain a complete observational census of baryons today? 

Here I have given a status report on these topics, which currently present more 
questions than answers. Major advances will only come through joint progress in 
observational capabilities, particularly in the ultraviolet, along with improving 
simulations for the evolution of baryons in the Universe. 
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